Universal Features of Information Spreading Efficiency on d-dimensional lattices 
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A model for information spreading in a population of -/V mobile agents is extended to d-dimensional 
regular lattices. This model, already studied on two-dimensional lattices, also takes into account 
the degeneration of information as it passes from one agent to the other. Here, we find that the 
structure of the underlying lattice strongly affects the time r at which the whole population has 
been reached by information. By comparing numerical simulations with mean-field calculations, 
we show that dimension d — 2 is marginal for this problem and mean-field calculations become 
exact for d > 2. Nevertheless, the striking nonmonotonic behavior exhibited by the final degree of 
information with respect to N and the lattice size L appears to be geometry independent. 
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I. INTRODUCTION 

The problem of information spreading among a popu- 
lation has been intensively studied in the lastyears and 
several aspects have been focused upon [H, H, [3, 13 • 

The population is generally represented by means of a 
graph such that an interaction (link) between two or more 
agents (nodes) means that there is a flow of information 
among them. Recently, the dynamics of agents making 
up the population has also been taken into account p], . 
Not only does the mobility of agents provide a realistic 
feature, but it also affects the network of acquaintances. 

In an earlier paper jl[ we introduced a model where 
agents are represented by N random walkers which dif- 
fuse on a square L x L lattice and possibly interact if 
they are close together. This model also takes into ac- 
count the degradation of information when passing from 
one agent to another: a decay constant z quantifies such 
alteration. As a consequence, the information spreading 
is a history-dependent process which is governed by the 
rules underlying the diffusion of N random walkers on 
the given space. 

Indeed, diffusion phenomena are dramatically affected 
by the topology of the underlying space [1, III IS @ • 
For example, on infinite lattices, the asymptotic proba- 
bility for a random walker to return to its starting point 
equals 1 in one and two dimensions, while for d > 3 there 
is a non-null probability to never return to the starting 
point. The former case is said to be recurrent, while the 
latter is called transient 0,0,0]. Moreover, in many pro- 
cesses concerning interacting random walkers, dimension 
2 plays the role of an upper critical dimension: it sepa- 
rates a higher-dimensional regime where the mean-field 
results are exact from a lower-dimensional regime where 
fluctuations become important. This is the case, for ex- 
ample, for two-species diffusion-limited reactions [![, or 
for the trapping of a random walker by diffusing traps 
@. In such processes, additional logarithmic corrections 
for the power laws in dimension d — 2 typically appear. 
Hence, in general, when dealing with diffusion one should 
also wonder how the laws describing the problem are af- 



fected by the geometry. 

The model introduced in Ref. [H for d = 2 is now 
extended to d-dimcnsional hypercubic lattices. Numer- 
ical simulations are carried for dimensions from d = 1 
up to d — 5. Analytical investigations are led which 
especially focus on the one-dimensional case and on a 
mean-field approach which provides good estimates for 
high-dimensional (d > 3) lattices. Most of the results 
presented here do also hold for general Euclidean (i.e., 
translationally invariant) lattices, since the large-scale 
topology of these systems (and quantities depending on 
it, such as the time r defined below in the low-density 
limit) depends solely on their dimension, and not on 
small-scale details. 

The main important quantities we are concerned with 
are the Population-Awareness Time r, which represents 
the average time necessary for the piece of information 
to reach the whole population, and the final degree of 
information per agent I ag (z). 

The time r depends on the system parameters N and 
L. Our numerical results show that in the low-density 
regime, and in every dimension d, this dependence can 
be factorized as t(N,L) = f(N)g(L), where / and g 
depend on d. Moreover, in the low-density regime, we 
find the asymptotic behaviors of both f(N) and g(L) to 
agree with mean-field calculations for d > 3, while for 
d = 2 deviations from the mean-field behavior appear 
and for d = 1 the results are radically different. We 
therefore argue that dimension d = 2 is marginal for 
the phenomenon under examination, and the mean field 
calculation of t is exact for d > 2. 

The most important result contained in Ref. [IJ con- 
cerns a nonasymptotic phenomenon: the nonmonotonic 
dependence of the final degree of information per agent 
T ag on N and L, with the emergence of extremal points. 
A process of optimization of the final information is 
therefore intrinsically nontrivial. We show here that 
the existence of extremal points is not a consequence of 
the special choice d = 2, but it arises in all dimensions 
d > 1 . It therefore appears as a universal and geometry- 
independent phenomenon, occurring at the crossover be- 
tween high- and low-density regimes. 
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The paper is organized as follows. Section |TT] is de- 
voted to the description of the model. Section IIIII con- 
tains analytical results; it is divided into high-density 
calculations (Sec. IIII A[) ; low-density calculations for 
d = 1 (Sec. IIII B[) : low-density calculations for d > 1 
(Sec. IIIICl) . Section ITVl shows results obtained by means 
of numerical simulations. We first consider the popula- 
tion awareness time r (Sec. IIV A[) . then the behavior of 
the final degree of information and the quantities that 
affects it fSecs. UVBl IP7U]) . Finally, Sec.|V]includes our 
conclusions and perspectives. 



II. THE MODEL 

The model analyzed in this work represents an exten- 
sion of the one introduced in an earlier paper [l[ . In this 
section we briefly recall how it works. 

We consider a population of TV random walkers 
(agents) moving on a <i-dimensional hypercubic lattice 
sized L and endowed with periodic boundary condi- 
tion. Agents are initially (t = 0) distributed randomly 
throughout the whole volume L d . We define the density 
of agents as p = N/L d ; the "low-density" regime is for 
p <C 1 and the "high-density" regime for p » 1. At 
each following instant each agent jumps randomly to one 
of the 2d nearest-neighbor sites. Notice that the same 
site can be occupied by more agents, i.e, there are no 
excluded- volume effects. 

At t — we assume that only one agent (called "Infor- 
mation Source" ) carries information, while the remaining 
N — 1 agents are unaware. Two agents can then inter- 
act if their distance on the underlying lattice is < 1 and 
if one of them is informed and the other unaware. By 
"interaction" we mean an information passing from the 
informed agent, say j, to the unaware one k with a fixed 
decay constant z (0 < z < 1): if j carries information Ij, 
then k becomes informed with information 1^ = z ■ Ii. 
Hence, the information carried by agent j is represented 
by the quantity Ij, < Ij < 1, and, in particular, when 
Ij > (Ij — 0) the agent is "aware" ("unaware"). 

Once an agent has become informed, it will never 
change nor lose its information. As a consequence, there 
exists a final time ta n at which the total information of 
the system can no longer evolve: at this time the infor- 
mation has reached every agent and the simulation stops. 
Such time is a stochastic quantity with average value r 
called the Population- Awareness Time (PAT) and stan- 
dard deviation denoted as a T . Part of this work is de- 
voted to studying the properties of r as a function of L 
and N. 

The total number of informed agents at a given time t 
is again a stochastic variable; we call n(t) its average over 
all the realizations of the system [n(0) = 1; n(oo) = N]. 
As a result of our model, the information carried by an 
agent is always a power z l of the decay constant, I being 
the number of passages from the Information Source to 
the agent. It is convenient to divide informed agents into 



levels, so that an agent belongs to level / when the in- 
formation it receives has undergone I passages from the 
Information Source and equals z l . We call n(l,t) the 
number of agents belonging to the Zth level at time t, av- 
eraged over all different realizations [n(t) = X)*=o ^)] - 
The average total information at time t is therefore the 
generating function of n(t), 

t 

z(m)=x>(m)« 1 . w 

2=0 

In particular, we are interested in the final degree of in- 
formation T(z), that is the total information achieved 
once the whole population has been informed, 

N 

l{z) = l(z, oo) = ^n{l,oo)z l . (2) 

1=0 

We also denote its average value per agent as I a g{z) = 
T(z)/N. The quantity n(l, oo) as a function of I is called 
the final distribution of the population on levels. 



III. ANALYTICAL RESULTS 

Although the model cannot be exactly solved in the 
general case, it is possible to provide approximate solu- 
tions in some limit cases. There are two time scales in- 
volved in the process: one for the diffusion of the random 
walkers and one for the information passing. When they 
are very different, approximate analytical approaches be- 
come feasible and give results in good agreement with 
the numerical simulations. When the two time scales 
are comparable, only a numerical approach is possible 
(Sec. E). 

In this Section we give analytical results for the PAT 
and the final distribution on levels in two limit cases. 
Section fill Al treats the high-density (p ^ 1) limit, with 
particular regard to the case d = 1. Section [ill Bl consid- 
ers the asymptotic low-density (p< 1) limit for d = 1. A 
general mean-field theory of this limit for all dimensions 
is given in Sec. IIII CI 

A. High-density regime 

When p > 1 , we can assume that the set of informed 
agents covers a connected volume of the lattice, and that 
this volume expands with a constant velocity (depend- 
ing on the density p and dimension d). We clarify this 
statement by considering d = 1. 

Let us consider a chain of finite length L, with N — > oo 
agents on it, and label the sites with the numbers from 
1 to L (Fig. [1]). The number of agents on a given site 
is p; for N — > oo (p — * oo), the probability that a given 
site is empty is 0. Also, for p — > oo, we assume that with 
probability 1 at least one of the p agents of a given site 
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FIG. 1: High-density approximation for d = 1. 
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FIG. 2: Low-density approximation for d = 1. 



will jump to the left and one to the right. Let the Source 
be in at t = 1: all the agents in 0, 1, and —1 will get 
informed. At t = 2 some of the newly informed agents 
jump on ±2; hence, the agents on ±2 and ±3 become 
informed. The motion of the information front decouples 
from the random motion of the agents; it expands with 
a deterministic law, with constant velocity of 4 sites per 
time unit. The time required to cover the whole chain is 
then 



r ~ L/4. 



(3) 



The border of the informed zone contains all the newly- 
informed agents, so each time step adds a new level and 
n(l, oo) = Ap for I > 1. This is the origin, for high 
densities, of the plateau observed in the one-dimensional 
distributions (Fig. [7]). 

For d > 1 it can be shown (see Ref. [lj for the case 
d = 2) that the volume of the informed zone is a d- 
dimensional polyhedron. The time it takes the border of 
the informed volume to reach the border of the lattice 
is again L/4, but now L/4 more instants are required to 
cover the rest of the lattice. Hence, in this case 



r ~ L/2. 



(4) 



The new agents added at each step cover a (d — 1)- 
dimensional surface, hence n{l, oo) ~ l d ^ 1 for I < L/4, 
and n(l - L/4, oo) ~ n(L/4, oo) - (I - L/A) d - X for I up 
to L/2. 



[2]). The average distance among two agents is p^ 1 = ^. 
Due to the low-density hypothesis, we can neglect inter- 
actions involving more than two agents. In this approach 
we divide the problem of diffusion among N agents into 
a sum of easier (three-bodies) problems. 

Let us consider the instant to when the rightmost agent 
ai informs an unaware agent a,2 (Fig. [2]). Let us call 03 
the next unaware agent on the right: the average distance 
from (Z3 to ai and 02 is L/N = p^ 1 (if p is small enough, 
we can consider a\ and a 2 to be on the same site). A 
calculation regardin g a n epidemic model in one dimen- 
sion similar to ours found that for low densities the 
velocity of the front propagation approaches p/2. Hence, 
the average time it takes one of the two aware agents to 
meet the unaware one is {L/N) /{p/2) = 2{L/N) 2 . If we 
now suppose that 03 has been first reached by 0,2, it will 
take again a time 2{L/N) 2 for one of them to reach the 
next unaware agent 0,4 on the right, and so on. There are 
about N/2 processes of this kind on the right-hand side 
and N/2 on the left-hand side, which provides 



N ' 



(5) 



Now, if a% belongs to level I (02 to level I + 1), 03 will 
belong to levels I + 1 or I + 2 with probabilities 1/2; 04 
will belong to levels / + 1 , / + 2 or I + 3 with probabilities 
1/4, 1/2, and 1/4, respectively, and so on. It is easy 
to show by induction, starting from a Source on level 
0, that at the ith information passing the new agent on 

i 



the left-hand side is on level I with probability 2 



(0 < I < i); the same for the new agent on the right-hand 
side. Hence, the average final number of agents on level 
I is 



N/2 

i{i,^)~Y, 2 ~ l 



(6) 



There is no easy closed form for this sum, but it can be 
plotted for any value of N: the curve displays a plateau 
of height 2, before decaying to 0. 

If we call AT{t) the increment of the total information 
at time t, we can write 



AI{i + 1) = Al{i 



z + 1 



and therefore 



X{z) 



z{z+\) 
2(1 -z) 



z + 1 



N/2 



(7) 



(8) 



B. Low-density regime in d = 1 

Let us consider a chain of length L and a population of 
N agents randomly distributed on it, with N <C L (Fig. 



C. Low-density regime in d > 1 

In the case of low density (p « 1) the time an in- 
formed agent walks before meeting an unaware agent be- 
comes very large. We adopt a mean-field approximation 
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by assuming that the agents between each event have the 
time to redistribute randomly on the lattice. In this ap- 
proximation, the probability that two given agents are 
in contact at a given time is pd = (t<i) , where (r^) is 
the average time for two random walkers to meet on a 
c?-dimensional cubic lattice. 

The process is an absorbing Markov chain with N states; 
the system is in state k when it has k informed agents. 
The chain starts from state 1 and evolves to the absorbing 
state (state N) . The transition matrix P can be written: 
the transition probability from a state A: to a state m as 
a function of N and pd is 



Pk 



N-k 

m — k 



I-(I-W)* 



m — k 



{i-PdY 



N-m 



for any N and pd- This is an upper triangular matrix, 

f N — k\ 

since the binomial coefficient , is for to < k. 

V m — k J 

We then make a low-density approximation: we expand 
matrix P to first order in pd to obtain 



Pk 



1 — m (N — to) \ 
m (N — m) pd 




for m — k, 
for m = k + 1 , 
elsewhere. 



This means that the system in the state to has a prob- 
ability 1 — m (JV — to) pd to stay in to and a probability 
to (N — to) pd to jump to state to+1. We now take matrix 
Q, the submatrix obtained from P subtracting the last 
row and column (those pertaining to the absorbing state), 
and compute the fundamental matrix F = (1 — Q) -1 ; a 
direct calculation shows that F is an upper triangular 
matrix given by 



where the Ud are dimension-dependent constants (for ex- 
ample, i*2 = 0.758...). As we will show in the follow- 
ing, the asymptotic dependence of r on L agrees with 
mean-field results for every d, while the breakdown of 
the mean-field theory shows up in the dependence on N 
for d < 2. 

It is possible to include the distribution on levels in the 
Markov chain analysis, but the calculations become very 
lengthy and we give only the final result. It is 



n(l, oo) 



1 



(JV-1)! 



|a(JVr,I + l)|, 



(11) 



and is exact for every N. Here, | • • ■ | denotes the abso- 
lute value, and s(to, k) is the Stirling number of the first 
kind (T3 |. The s(m, k) are integers that appear in many 
combinatorial problems; one of the possible asymptotic 
expansions for Stirling numbers is [131 ] 



so that 



n(l, oo) 



ln(N) 1 
11 



(12) 



which is the form we will use to fit the low-density dis- 
tributions. From this distribution it also follows that 



1 - N z 



IV. NUMERICAL RESULTS 



A. Population- Awareness Time 



(13) 



I — rrr- — \ — > for A: > to 

Fu = < m (-™- m )Pd' — 

1 for k < to. 

The mean time r required to reach the absorbing state 
N starting from state 1 is given by the sum of the first 
row of F, 



1 

Pd 



N-l 

E 



1 



' m(N — m) 



(9) 



and for N — > oo 
2 



Np d 



( 7 + ln(jV))=2(r rf ) 7+ y jV) , (10) 



where 7 = 0.577... is the Euler-Mascheroni constant. 

A classical result [1, [l(| states that for d-dimensional 
cubic lattices the asymptotic dependence of {to) on the 
lattice size L is 



iti L 2 d=l 
u 2 L 2 ln(L) d = 2 
u d L d d>2, 



In this section we focus on numerical results concerning 
the Population- Awareness Time r. We recall that r has 
been defined as the average time it takes the piece of 
information to reach the whole population. Due to the 
analytical results discussed in the preceding section, we 
expect the functional form displayed by r(N, L) to be 
strongly affected by the topology of the lattice underlying 
the propagation. 

Figures [3] and [4] show the dependence of r on L with 
N fixed, and on N with L fixed, respectively. In Fig. 
where both the results for d — 1 and d = 3 are displayed, 
two different regimes, of low and high density, are clearly 
distinguishable for dimension d = 1; for d = 3 the range 
of the high-density regime is too small, and only the low- 
density one can be seen. 

The high-density behavior is independent of N; indeed, 
for p ^> 1 we find 



T=-,d=l 

4' 



T = -,d>l, 



(14) 



in agreement with Eqs. (O and ([!]). 
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In the low-density regime (p <C 1) and for d ^ 2, r 
follows the behavior calculated in Sees. IIII "Bl and [Til CI 
For d = 1 the results as a function of L and N are fitted 
by 



dN a L p , 



(15) 



with Ci = 0.96(5); a = -0.98(5); (3 = 1.98(2), in agree- 
ment with Eq.©. 

For d — 3, 4 we found the best fit for r to be given by 
the function 



C d 



ln(L) +A 
N 



(16) 



where (3 = d within a 1% error and A = 0.59(3), in agree- 
ment with the value 7 = 0.577 . . . found in Eq. (JTOJl . The 
values of the constants are C3 = 0.77(1); C4 = 0.39(1) 
(different in general from the Ud of the mean-field ap- 
proximation) . 

The low-density limit in the case d — 2 deserves a 
separate discussion. It is still possible to express r as a 
product of two distinct functions, 



t - f 2 (N) L 2 ln(L), for d = 2. 



(17) 



The dependence on L is in agreement with the improved 
mean- field calculation (see Sec. IIII C[) . This best fit is 
better than that in Ref. [l[, where we hypothesized a 
non-integer power law (L 2 2 ). 

The analytical form of f2(N) cannot be unequivocally 
determined by the simulations. In the fitting range the 
function ^M-^) agrees with the numerical results better 
than the power-law TV" - 66 previously given However, 




FIG. 3: Log-log scale plot of the Population- Awareness Time 
r versus the lattice size L. Results obtained for a chain 
(dashed line) and for a cube (solid line) are depicted. Dif- 
ferent lattice-size values are shown with different symbols, as 
explained by the legend. For large densities (p S> 1, HD) and 
low densities (p <C 1, LD), straight lines represent the best fit 
according to Eq. (|14[) and Eqs. ( I15|) and (|16[) . respectively. 
Error on data points is less than 2%. The standard deviation 
is not appreciable on this scale. 



in this case the value of the fitting parameter is A — 
—0.18(2), hence is definitely different from the mean-field 
value 7. Since there are no analytical calculations to 
support this functional form with this particular value of 
the fitting constant for d — 2, we cannot rule out higher- 
order logarithmic corrections. 

To summarize, in the low-density regime, the func- 
tion t(N, L) factorizes into two parts, depending, respec- 
tively, on L and N: 



L 2 



d= 1, 



/ 2 (iV)L 2 m(L), d = 2, 



(18) 



a 



l + HN) Td 



N 



■ L 



the Cd being dimension-depending constants. The most 
satisfying fitting function we have found for fz(N) is 



A+ln(N) 
N 



A 



-0.18. 



The standard deviation <r T displays a similar depen- 
dence on N and L for low densities: <j t cc N^L 2 for 
d = 1, and so on. For high densities o> becomes van- 
ishingly small, which is explained by the fact that the 
propagation of information becomes a deterministic pro- 
cess. 



Final Distribution on Levels: Universality of 
the Extremal Distribution 



In Sec. |TT] we introduced the function n(l, 00), which 
represents the final distribution of agents on levels and 
is strongly connected with the final degree of informa- 
tion T{z). The asymmetrical-bell shape displayed by 
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FIG. 4: (Color online) Dependence of the Population- 
Awareness Time r on the number of agents TV for the cubic 
lattice d = 3; different values of the lattice size L are shown 
with different symbols and colors. When the density of the 
system is low, data points lie on the curve given by Eq. (|16p 
which represents the best fit. Error on data points is less than 
2%. 
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FIG. 5: (Color online) Final population distribution on levels 
nil, oo) for low-density systems. Data points agree with the 
fitting line drawn according to Eq. (|19[) . For d = 2 the 
fitting parameters depend smoothly on both L and N, and 
the curve is distinct from those of higher dimension and same 
AT. For d > 3, systems of different dimension d and size L 
display distributions that overlap within the error. Only the 
dependence on N is left, as is shown for d = 3, iV = 512. 
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FIG. 6: Final population distribution on levels for high den- 
sities and d = 1,2,3. The dependence on I is a power law: 
n(l, oo) ~ Z d_1 . 



the distributions for hypercubic lattices with dimension 
d > 2 (Fig. [5]) and the way they evolve while varying 
the system parameters N and L are analogous to the 
two-dimensional case [l[. 

In the limit case p> 1 (Fig. [6} and in every dimension, 
the final distribution on levels follows the law n(l, oo) ~ 
Z d_1 , in agreement with the calculation in Sec. IIII A\ 

For p small enough (i.e., for L > L and N < N, see 
below) the population distribution on levels for d > 2 is 
well fitted by the following function: 



n(l, oo) 
N 



A 



(In Ny 
T(Bl + C)' 



(19) 



where T(x) is the Euler gamma function (Fig. [5]). The 
previous equation is a generalization of Eq. (|12|) found in 
the mean-field approximation for the low-density regime. 

The fitting parameters A,B,C at low densities for d = 
2 are smoothly dependent on N and L, while those for 
d > 2 are all close to 1 and independent of the lattice 
size L, keeping only the dependence on N. Moreover, 
the data points for d > 3 with the same values of N all 
collapse on the same curve. This holds for all the regular 
lattices with d > 3 we have considered: the distribution 
curves at low densities are independent of d, and agree 
with the mean-field form [Eq. (|12| . 

The description given so far concerns d > 2 lattices; 
in Fig. [7| we show for N = 512 and varying L the one- 
dimensional case, which exhibits quite different distribu- 
tions. Such distributions are still very sharp for very high 
values of the density p, but soon develop a plateau by in- 
creasing L; the plateau persists up to low densities. The 
existence of a plateau was justified both in a high-density 
and in a low-density approximation (Sec. IHIjl . 

The main result of Ref. 1] concerned the existence of 
an extremal curve for the distribution of agents on lev- 
els. We have found that this feature does not depend on 
the dimension d of the lattice. We show in Fig. [7] how 
the extremal distribution emerges in d = 1, as a func- 
tion of L, for a particular value L = L (here, L ~ 1024), 
and keeping N fixed, notwithstanding the fact that its 
shape is dramatically different with respect to the higher- 
dimensional ones. While L < L, the distribution dis- 
plays a plateau whose height (width) is a monotonically 
decreasing (increasing) function of the chain length L; 
the distribution curve shifts to the right. Conversely, for 
L > L a shift-back phenomenon analogous to that dis- 
cussed in Ref. [l[: now, by rising L, the height gets larger 
while the width gets smaller. As can be seen from Fig. 
[31 L corresponds to the crossover between high- and low- 
density regimes. The same happens by varying N and 
keeping L fixed; there is an extremal distribution for a 
particular value N, depending on L, and corresponding 
to the crossover between the two regimes. 

This shift-back phenomenon, and the existence of an 
extremal distribution, occur in all the dimensions we have 
investigated (up to d = 5). It therefore constitutes a 
universal feature, independent of lattice dimension, and, 
as we will see, it provides striking effects on the final 
degree of information. 



C. Degree of Information 

In this section we deal with the final degree of infor- 
mation T(z) — I{z, oo) ([2]) and its dependence on the 
decay constant z and system parameters N, L. We re- 
mind [Eq. (JTJ>] that T(z) is the generating function of the 
final populations n{l, oo), hence its value depends on the 
final distribution of the population on levels analyzed in 
the preceding paragraphs. 

Let us firstly consider the dependence on the decay 
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FIG. 7: (Color online) Population distribution on levels at 
t — t for one-dimensional systems with N = 512 and L rang- 
ing from 2 3 to 2 14 , as shown by the legend (the lines are guides 
to the eye). The behavior of the distribution is nonmono- 
tonic with respect to L: by increasing L from small values, 
the curves first shift to the right and flatten; the rightmost, 
extremal curve corresponds to L — 1024. Then, by further in- 
creasing L, the curves shift back to the left and sharpen. The 
inset shows in detail the shift back with the curves pertaining 
toL = 2 8 ,L = 2 10 ,L = 2 14 . 

constant z. Again, results highlight strong differences be- 
tween the one-dimensional and higher-dimensional cases 
(d > 2). In the latter case and for the low-density regime 
(approximately p < 2~ 8 ), we find 

l{z) = N z , (20) 

within the error (< 4%). 

On the other hand, when d = 1, the final degree of 
information shows an exponential growth which can be 
represented by the following equation: 

I(z)^A Z ^^\ l-e- B - N ^), (21) 

where A and B smoothly depend on N and L. Equa- 
tions (|20"|) and (|2T|) are in very good agreement with the 
expressions found in the mean- field approximation [Eqs. 
(fT^I and ©, respectively]. 

Once z is fixed, X ag (z) depends nonmonotonically on 
N and L: let us follow it for N fixed and varying L in 
Fig. [8] in the two cases d = 1 and d = 3. For L small, 
due to the narrow distribution discussed in the preced- 
ing section, the value of the information is high. When 
L = L, the population distribution on levels reaches its 
extremal form and the information displays a minimum. 
As L increases, the information starts to rise again (as 
can be seen, the effect gets more marked by increasing 
the dimension). Hence, given a population number N, 
there is an optimal lattice size L for which the final in- 
formation is minimum. The same happens having fixed 



L and letting N vary: there is a minimum for N = N, 
depending on L. As underlined in Ref. [l[, the existence 
of a local minimum of the final information implies that 
choosing an optimization strategy for the spreading of 
information on the lattice is not trivial. There is no a 
priori right direction in parameter space where to move 
in order to improve I(z); rather, the direction depends 
on the starting point. 

V. CONCLUSIONS AND PERSPECTIVES 

In this work the model of information spreading previ- 
ously introduced has been extended to different geome- 
tries; indeed, we considered the chain and d-dimensional 
hypercubic lattices. The occurrence of a nonmonotonic 
behavior for the final degree of information is not due to 
a special geometry underlying the process, but its origin 
lies in the crossover between the two different regimes of 
high and low density. Therefore, the existence of min- 
ima in the final degree of information is universal and, 
remarkably, even the possibility to derive optimization 
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FIG. 8: (Color online) Semilog scale plot of final degree of 
information per agent 1 ag (z) = I(z)/N vs lattice size L, for 
d — 1 (top) and d — 3 (bottom). The decay constant is fixed 
at z = 0.95. Several values of N are shown with different 
symbols and colors (lines are guides to the eye) and the legend 
is the same for both figures. Notice that minimum depth is 
greater for the latter case. Error on data points is < 4%. 
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strategies does not depend on the particular structure 
the process is embedded in. 

On the other hand, the asymptotic laws for r are inter- 
estingly related to the geometry underlying the random- 
walk diffusion. In particular, d — 2 is a marginal dimen- 
sion separating two well-behaved cases, which suggests 



an investigation on in-between dimensions [14j . 

The robustness of the existence of extremal point for 
X is an important point since the possibility of extract- 
ing optimal strategies is not a feature restricted to some 
special structures. 



[1] E. Agliari, R. Burioni, D. Cassi and F.M. Neri, Phys. 

Rev. E 73, 46138 (2006). 
[2] M.C. Gonzalez, P.G. Lind and H.J. Herrmann, Phys. 

Rev. Lett. 96, 88702 (2006). 
[3] L. Huang, K. Park and Y. C. Lai, Phys. Rev. E 73, 35103 

(2006). 

[4] R.V. Sole and S. Valverde, Physica. A 289, 595 (2001). 
[5] E.W. Montroll and G.H. Weiss, J. Math. Phys 6, 167 
(1965). 

[6] R. Burioni and D. Cassi, J. Phys. A 38, R45-R78 (2005). 
[7] G. Polya, Math. Ann. 84, 149 (1921). 
[8] U. C. Timber, M. Howard and B. P. Vollmayr-Lee, J. 
Phys. A 38, R79-R131 (2005) (1965). 



[9] F. van Wijland, J. Phys. A: Math. Gen. 35, 5391-5401 
(2002). 

[10] Indeed, the result in [f| concerns the mean absorption 
time into a trap at the origin for a single walker. How- 
ever, the problem of the mean first meeting time of two 
walkers, 1 and 2, can be mapped onto the problem of 
absorption at the origin with a proper rescaling of time. 

[11] C. P. Warren, G. Mikus, E. Somfai, and L. M. Sander, 

Phys. Rev. E 63, 056103 (2001). 

[12] http://mathworld.wolfram.com/StirlingNumberoftheFirstKind.html 

[13] H. S. Wilf, J. Comb. Th. Ser. A, 64, 344 (1993). 

[14] E. Agliari, R. Burioni, D. Cassi and F.M. Neri, in prepa- 
ration. 



